Abstract

Introduction

Knowing species interactions

Knowing who eats whom in a a food web can help us answer crucial ecological questions such as: How would the food web respond if the top predators are removed? Which species is a keystone species in the food web? Which species should be given priority when conservation efforts are implemented? (citreqd)

Conventional methods of knowing interactions

These species interactions can be inferred from different types of information and methods such as gut content analysis (Peralta-Maraver, L’opez-Rodr’ıguez, and de Figueroa 2017), stable isotope ratio composition of tissues (Layman et al. 2007), experimentation (Warren 1989), DNA metabarcoding (Roslin and Majaneva 2016). These methods vary in the type and amount of food web information that can be extracted. For e.g.: gut content data of a species reveals which prey species are fed by that predator species in the food web while stable isotope ratio data of a species tells us the approximate trophic position of that species in the food web. As a food web can be a complex network with multiple links, it requires a lot of data which means intensive sampling. It might not be feasible to use a single method to extract full information about the food web interactions.

  • Food web properties are highly sensitive to varying levels of sampling effort.

  • How do we quantify the sampling effort?

Gut content analysis

Gut content analysis is performed to know what an individual has eaten, and the items in the stomach are used to infer the diet of that species. There is considerable uncertainty involved in this method (Baker, Buckland, and Sheaves 2014). In fish, we find ambiguity while separating prey items in the gut. There are loose tissues which are not identifiable and cannot be assigned to a specific prey item with certainty. Even when it is possible to separate prey items and identify the prey, there are lot of factors such as sample size of consumers, mechanical prey handling, differential digestion and evacuation rates of different types of different prey types and volumes, and the ingestion order that in combination result in unquantifiable error which are difficult to interpret in the predator diet (Hyslop 1980; Rindorf and Lewy 2004).

Stable isotope ratios

Stable isotope ratios (SIR) of tissues from a single individual is used to quantify diet over varying temporal and spatial scales from one sample (Crawford, Mcdonald, and Bearhop 2008). Using SIR, several different sources of information can be extracted from a single sampling event as different tissues integrate diet over different temporal scales as well as spatial scales. However, there is uncertainty involved in this method which is due to variability in the isotopic fractionation values across multiple combinations of diets and tissues/species, unquantified temporal or spatial variation in prey isotopic values and variation caused by routing of particular dietary nutrients into particular tissues.

Sampling effort problem

These methods of inferring prey predator interactions require a lot of sampling effort (citreqd). For e.g.: one needs to know the information of the gut of all the species to construct the food web using only gut content data. The amount of sampling effort affects how well the food web properties are predicted (Goldwasser and Roughgarden 1997). The amount and type of information provided by these methods vary and some of the information inferred might overlap among these methods. The question we would like to address here is how should the sampling effort be distributed among the methods/species? Can we infer the complete structure of a food web using information from only one method? For e.g.: Is it possible to know all the prey predator links by using only the stable isotope ratio data? Do we need to sample all the species?

General research questions

  • How much gut content data is required to know the food web structure?
  • How do we combine different types of food web data for food web prediction?
  • How do we incorporate uncertainty in food web data to food web prediction?

Research questions specific to gut content data

  • Can ADBM predict the food web structure using partial gut content data?
    • first, we use a simulated food web and simulated gut content data to answer it.
    • second, we use real gut content data from Broadstone Stream data to answer it
  • How does the result from above questions generalise over food webs?

Our approach

An important question in food web ecology is to ask how much data of prey-predator interactions is required to know the food web structure (citreqd). Is it possible to construct the food web by knowing only partial data? How should the sampling effort be distributed across species? Here, we present a theoretical study as well use a real data to answers these questions. We use a mechanistic model known as allometric diet breadth model (ADBM) which predicts food web structure. We provide partial data to the the model and observe if we can predict the food web structure and also quantify the accuracy of the prediction. We show how one can use multiple sources of food web data and infer the structure of a food web.

Using real gut content data and stable isotope ratios data we predict the Broadstone Stream food web.

Why gut content data and SIR data only?

  • their importance
  • amount of information
  • ease of extraction
  • availability of data

Why ABC?

To address the above questions, we use a Bayesian approach known as approximate Bayesian computation (ABC). With the help of summary statistics in ABC, we can use the partial and diverse types of information to predict the structure of the food web.

Materials and Methods

In the upcoming sections, we present a detailed description of how we synthesize gut

data synthesis and ADBM’s predictions

First, we address these questions using in silico (Fig ) by simulating food web structure from the allometric diet breadth model (ADBM) whose true parameters are known. We then synthesize data such as gut contents, stable isotope ratios, direct observation from the simulated food web. Next step would be to incorporate uncertainty in these data so that it mimics similar challenges as in real data. Then, we use these observed data which includes the uncertainty to predict the food web from the ADBM parameterised using approximate Bayesian computation (ABC). Here, we make use of the summary statistics in the rejection ABC approach.

\label{fig:fig_m1}  Flowchart showing the method

Flowchart showing the method

Gut contents simulation

Stable Isotope Ratios simulation

Broadstone Stream example

  • Broadstone Stream (51deg 05min N 0deg 03min E; 120 m above sea-level) is a second-order tributary of the River Medway in south-east England (see Hildrew (2009) for a detailed site description)
  • The acidity of the stream (pH 4.7–6.6) excludes fish, resulting in an invertebrate-dominated food web (Woodward and Hildrew 2001). There are about 25 common invertebrate species (Woodward, Speirs, and Hildrew 2005). Among the common predators, there are three large species (Cordule- gaster boltonii Donovan, Sialis fuliginosa Pict. and Plectrocnemia conspersa [Curtis]) and three small species (the larvae of the tanypod midges Macrope- lopia nebulosa [Meigen], Trissopelopia longimana [Staeger], and Zavrelimyia barbatipes [Kieffer]), which we focus on in detail here when investigating size- structuring within the intraguild subweb.

ABC

Result

Prediction using simulated gut content data from a simulated food web

\label{fig:fig_r0} True skill statistics between predicted food web and observed food web for simulated small reef food web estimated for different number of predator species' guts sampled. For every value in the horizontal axis, there are 100 samples of randomly chosen species from a set of 50 species whose gut content data are simulated.

True skill statistics between predicted food web and observed food web for simulated small reef food web estimated for different number of predator species’ guts sampled. For every value in the horizontal axis, there are 100 samples of randomly chosen species from a set of 50 species whose gut content data are simulated.

We investigated how the true skill statistics (TSS) between the predicted food web and observed food web varied with the number of predator species’ sampled (Fig ) for gut data. This investigation was done for simulated small reef food web using the ADBM. For every n number of predator species, we had 100 sets of gut data randomly selected from 50 species, where a set consists of n number of different predator guts. The TSS quickly saturated with increasing number of predator species. When compared with the partial individual gut data, the TSS estimated from that of the complete individual gut data had less uncertainty and higher mean TSS. The maximum values of TSS for complete gut data and partial gut data were almost equal, with the minimum TSS values for partial gut data lower than that of the complete gut data. Eventually, the gap between the mean TSS for partial gut data and complete gut data reduced with increasing number of predator species suggesting when we have enough predator species’ gut data, the achieved TSS was almost constant whether the gut data of an individual predator was complete or not. The maximum TSS for complete gut data was very close to one and almost remained constant wih increasing number of different predator species sampled. With gut data of only five predator species, 95% of the maximum mean TSS was achieved for complete gut data, while with partial gut data the maximum mean TSS was achieved with 15 predator species. This shows that one does not need to know the gut data of all the species to predict the food web and the number of predator species is dependent on the completeness of an individual gut data. Hence, it is not necessary to know 100% accurate information of the gut data.

Prediction using empirical (real) Broadstone stream gut content data

\label{fig:fig_r1} True skill statistics between the predicted and observed Broadstone stream food web computed for number of different predator guts sampled (real gut data and randomly generated gut data). There were a total of 1008 predator guts. Blue curve corresponds to real gut data and red curve corresponds to randomly generated gut data. Line represents the mean and shaded region represents the range/uncertainty. Every point on the horizontal axis correspond to 100 samples.

True skill statistics between the predicted and observed Broadstone stream food web computed for number of different predator guts sampled (real gut data and randomly generated gut data). There were a total of 1008 predator guts. Blue curve corresponds to real gut data and red curve corresponds to randomly generated gut data. Line represents the mean and shaded region represents the range/uncertainty. Every point on the horizontal axis correspond to 100 samples.

Prediction using simulated gut content data from a simulated food web

\label{fig:fig_r2} True skill statistics between predicted food web and observed food web for simulated small reef food web estimated for different number of predator species' guts sampled. For every value in the horizontal axis, there are 100 samples of randomly chosen species from a set of 50 species whose gut content data are simulated.

True skill statistics between predicted food web and observed food web for simulated small reef food web estimated for different number of predator species’ guts sampled. For every value in the horizontal axis, there are 100 samples of randomly chosen species from a set of 50 species whose gut content data are simulated.

The value of true skill statistics between the predicted food web and observed food web quickly saturated with increasing number of sampled real predator guts (Fig ). With only 228 number of real predator guts, 90% of the maximum mean TSS was achieved while 95% of the maximum TSS was achieved with around 428 real gut content data. The uncertainty (shaded region) in the TSS reduces with increasing number of predator guts. Although at low number of predator guts the maximum TSS is quite high but the uncertainty is high as well. This is quite possible just by chance which is suggested by the overlap of TSS region with that of randomly generated gut content data.

We compared the TSS predicted from real gut content data with randomly generated gut content data. Unsurprisingly, the range/uncertainty of TSS was high and almost remained constant with increasing number of predator guts sampled. The mean TSS fluctuated for lower number of predator guts sampled but eventually became constant, and was close to 0 which suggests the prediction is no better than random.

We investigated how the TSS varied with the number of predator guts sampled in fig where the predator guts were divided into two cases depending on whether that predator has body size less than or greater than the average body size. We found a saturating curve in both the cases. The TSS corresponding to gut data less than average body size was higher than that corresponding to gut data greater than average body size for any number of predator guts sampled. Also, the TSS corresponding to the gut data less than average body size was higher when compared with the TSS corresponding to the gut data with all the body sizes.

\label{fig:fig_r3} True skill statistics between the predicted food web and observed food web estimated with respect to the number of predator guts. The estimation is done for three sets of gut data: gut content data of predators less than the global mean body size, greater than the global mean body size and all the gut content data.

True skill statistics between the predicted food web and observed food web estimated with respect to the number of predator guts. The estimation is done for three sets of gut data: gut content data of predators less than the global mean body size, greater than the global mean body size and all the gut content data.

Prediction using simulated stable isotope ratio data from a simulated food web

\label{fig:fig_r4} True skill statistics between predicted food web and observed food web for simulated small reef food web estimated for different number of predator species. For every value in the horizontal axis, there are 100 samples of randomly chosen species from a set of 50 species whose gut content data are simulated.

True skill statistics between predicted food web and observed food web for simulated small reef food web estimated for different number of predator species. For every value in the horizontal axis, there are 100 samples of randomly chosen species from a set of 50 species whose gut content data are simulated.

Combining both gut content data and stable isotope raios data

Discussion

Include digestion in gut content simulation

extrapolating ADBM’s predictions to general food web

We first answer the questions with the help of ADBM. Now, the question is how much does it translates to real food web.

extension to other food web models

Conclusion

Appendix

References

Baker, Ronald, Amanda Buckland, and Marcus Sheaves. 2014. “Fish Gut Content Analysis: Robust Measures of Diet Composition.” Fish and Fisheries 15 (1): 170–77. https://doi.org/10.1111/faf.12026.

Crawford, Kerry, Robbie A. Mcdonald, and Stuart Bearhop. 2008. “Applications of Stable Isotope Techniques to the Ecology of Mammals.” Mammal Review 38 (1): 87–107. https://doi.org/10.1111/j.1365-2907.2008.00120.x.

Goldwasser, Lloyd, and Jonathan Roughgarden. 1997. “Sampling Effects and the Estimation of Food-Web Properties.” Ecology 78 (1): 41–54. https://doi.org/10.1890/0012-9658(1997)078[0041:SEATEO]2.0.CO;2.

Hildrew, Alan G. 2009. “Chapter 4 Sustained Research on Stream Communities.” In Advances in Ecological Research, 41:175–312. Elsevier. https://doi.org/10.1016/S0065-2504(09)00404-8.

Hyslop, E. J. 1980. “Stomach Contents Analysis—a Review of Methods and Their Application.” Journal of Fish Biology 17 (4): 411–29. https://doi.org/10.1111/j.1095-8649.1980.tb02775.x.

Layman, Craig A., D. Albrey Arrington, Carmen G. Montaña, and David M. Post. 2007. “Can Stable Isotope Ratios Provide for Community-Wide Measures of Trophic Structure?” Ecology 88 (1): 42–48. https://doi.org/10.1890/0012-9658(2007)88[42:CSIRPF]2.0.CO;2.

Peralta-Maraver, I., M. J. L’opez-Rodr’ıguez, and J. M. Tierno de Figueroa. 2017. “Structure, Dynamics and Stability of a Mediterranean River Food Web.” Marine and Freshwater Research 68 (3). CSIRO PUBLISHING: 484–95. https://doi.org/10.1071/MF15154.

Rindorf, A, and P Lewy. 2004. “Bias in Estimating Food Consumption of Fish by Stomach-Content Analysis” 61: 12.

Roslin, Tomas, and Sanna Majaneva. 2016. “The Use of DNA Barcodes in Food Web Construction—Terrestrial and Aquatic Ecologists Unite!” Edited by Elizabeth Clare. Genome 59 (9): 603–28. https://doi.org/10.1139/gen-2015-0229.

Warren, Philip H. 1989. “Spatial and Temporal Variation in the Structure of a Freshwater Food Web.” Oikos 55 (3): 299. https://doi.org/10.2307/3565588.

Woodward, Guy, and Alan G. Hildrew. 2001. “Invasion of a Stream Food Web by a New Top Predator.” Journal of Animal Ecology 70 (2): 273–88. https://doi.org/10.1111/j.1365-2656.2001.00497.x.

Woodward, Guy, Dougie C. Speirs, and Alan G. Hildrew. 2005. “Quantification and Resolution of a Complex, Size-Structured Food Web.” In Advances in Ecological Research, 36:85–135. Elsevier. https://doi.org/10.1016/S0065-2504(05)36002-8.